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We consider here the problem of a "giant spin", with spin quantum number 5* 2> 1, interacting 
with a set of microscopic spins. Interactions between the microscopic spins are ignored. This model 
describes the low-energy properties of magnetic grains or magnetic macromolecules (ferromagnet- 
ically or antiferromagnetically ordered) interacting with a surrounding spin environment, such as 
nuclear spins. 

Our aim is to give a general method for truncating the model to another one, valid at low 
, energies, in which a two-level system interacts with the environmental spins, and higher energy 

' terms are absorbed into a new set of couplings. This is done using an instanton technique. 

We then study the accuracy of this technique, by comparing the results for the low energy effective 
■ Hamiltonian, with results derived for the original giant spin , coupled to a macroscopic spin, using 

exact diagonalisation techniques. We find that the low energy central spin effective Hamiltonian gives 
very accurate results (with increasing accuracy for large S), provided the typical coupling energies 
^ ■ between the giant spin and the microscopic spins are not too large, and provided temperature and 

0^ ' external field are sufficiently low. The essential limitation to the applicability of the low-energy 

effective Hamiltonian is just the semiclassical WKB approximation itself, which inevitably fails for 
very small 5*. 

Our results thus justify previous use of this effective Hamiltonian in calculations of the effects of 
nuclear spins on the dynamics of nanomagnetic systems. 



00 

o 



B 

o 



I. INTRODUCTION 



In this paper we examine a model which was originally introduced to deal with the effects of environmental spins 
. on the coherent tunneling of magnetic grain magnetisation It consists of a central "giant spin" having 

I spin quantum number 5 1, coupled to a set {(7k} of "microscopic spins"; the microscopic spins are only very 
weakly coupled to each other. The model is then completely described by a Hamiltonian Ho{S) for the central spin, 
,({(7fc}) for the microscopic spins, and a coupling Hint{S, {<?*:}) between the two: 



1 ^ 

H{S,au) ^ Ho{S) + — Y,>^kS ■ ak + Vd,j>{{ak,ak'}) , (1.1) 
fe=i 

in which we specialise immediately to the form relevant to most of this paper, in which the interaction between the 
central spin and the environment has a hyperfine form, and the dynamics of the environmental spins (considered 
now to be nuclear spins), comes from their mutual dipolar or indirect Suhl-Nakamura interactions. In most practical 
applications of this "central spin model", the environmental nuclear spins will be either inside the object carrying the 
giant spin, or near it, in some substrate, or solvent, or surrounding matrix. The central spin may be a magnetic grain 
or a magnetic macromolecule such as ferritin |Q or, on a smaller scale Mni2-ac. Similar models were also 
introduced some time ago to describe the large superparamagnetic " spin clusters" which are believed to exist in many 
disordered magnets at low temperature, such as Si : P near the metal- insulator transition ||^, or "giant magnetic 
polarons" [ p^ . Similar spin clusters exist in " quantum spin glasses" pT| . In these models the giant spin is somewhat 
less uniquely defined (indeed its size depends on temperature). Finally, a similar central spin model was considered 
a long time ago by Gaudin, who wrote down a set of Bethe ansatz equations for it p^ . 
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We should also note that models like (1.1) can be of interest even for non-magnetic systems (although the coupling 
term will not then be a hyperfine one). This is because most large quantum objects will have (at least) a coupling 
to their own nuclei, and perhaps to other quasi-free spin degrees of freedom like paramagnetic impurities. Thus in 
the oft-cited example of the superconducting SQUID ring , there is an electromagnetic coupling of the tunneling 
flux coordinate (and the associated ring supercurrent) to all of the nuclear spins within a penetration depth of the 
surface (either in the junction or in the ring), as well as to any paramagnetic impurity spins that happen to be 
present jl^. The nuclear spins can then not only change the energetics of the relevant collective coordinate, but also 
dephase it, in a kind of "inverse Stern-Gerlach effect", where nuclear spins play the role of measuring devices | p^ , p^ . 
It then becomes important to determine how the coupling to the background spins will affect the quantum dynamics 
of the central system. This depends of course on the system; for example, in the case of SQUID's, it is fairly easy to 
show that the effect of the environmental nuclear and paramagnetic spins on flux tunneling is usually negligible, but 
that the effect on "Macroscopic Quantum Coherence" may not be. On the other hand, the strength of the hyperfine 
coupling, between the order parameter of many nanomagnetic systems and their nuclear spins, is such that the nuclear 
environment can seriously disrupt the quantum dynamics of the nanomagnet at low T . 

The purpose of the present paper is not to analyze the effects of the environmental spin bath on the dynamics of 
nanomagnet - we have done this in a number of other papers. Instead we wish to investigate the validity of the effective 
Hamiltonian we have previously employed. This Hamiltonian results from a truncation to a lower-dimensional Hilbert 
space in which the giant spin dynamics are reduced to that of a two-level system (still coupled to N environmental 
spins). This truncation is analogous to the truncation of a quantum system like a SQUID, coupled to delocalized 
electrons, to the "spin-boson" model at low T. In the present case of the central spin, the truncated central spin 
effective Hamiltonian takes the form 



iJef f ~ \ 2Aof_ cos 



7r5 -f ^ aun ■ Uk ~ Pon ■ Ho 



H.c. 



, N ^ N N N 

y ^ w| ffe • O'fc + - ^ TOfc • <?fc + ^ ^ K'fc'? ^k^k' 

k=l k=l k=lk' = l 



(1.2) 



In this Hamilto nian the two-level dynamics of S will be described by the Pauli vector matrix r , and we have assumed 
that Ho{S) in ( |l.l| ) contains an external field Ho- This is the same form as we have used in previous investigations 
1^,^. Here Aq is a tunneling matrix element, n, Ik, and ifik are unit vectors, with and rhk mutually perpendicular; 
Q!fc and Po are dimensionless parameters, and wjl and lo^ are "longitudinal" and "transverse" diagonal coupling energies 
0- 

In the present pap er we will show how one goes from (LI) to ( [l.2| ), in an instanton calculation; and we w ill d iscuss 
how accurately (L^) represents the dyn amic s of (IT), by carrying out exact diagonalization studies of (IT), and 
comparing them with the predictions of (O). To do this we will pick a special (but actually quite generic) form for 
Ho{S), viz., a biaxial (easy axis/easy plane) form 



Ho{S) 



Kn Si + K^Sl 



(1.3) 



and calculate the parameters in (1.2) deriving from ( p_.3| ), in an instanton expansion. We should note here that some 
similar exact diagonalisation studies of this model were recently carried out by Abarenkova & Angles d'Auriac [pT| , 
who also compared their exact diagonalisation results with a somewhat simplified version of (1.2) in certain limits. 
In the present work we have a rather different aim; we wish to show that the structure of the exact diagonalisation 
results can be put in a form very close to that in ( |l.2D , and then we can check numerically whether the coefficients 
in ( |l.2D agree with those we calculate using an instanton procedure, over a very wide paramenter range. Thus we are 
not just trying to test the truncated model numerically; we are also trying to test its form, and to see if anything is 
missing from it. 



We find that, as assumed in previous work, 1.2) gives a very good description of the low-energy dynamics of 
(1.1), both as regards the mathematical form and the numerical values calculated for the coefficients in an instanton 
expansion, provided that 

(i) The semiclassical approximation for the "bare" giant spin Hamiltonian Ho{S) is a good one; 

(ii) The parameter Eo — uJqN^^'^, characterizing the spread in energies caused by the hyperfine coupling, and 
the thermal energy ksT, satisfy the inequality Eo, ksT <C ^^07 where VLo is the "bounce frequency" associated with 
semiclassical tunneling, which characterizes roughly the energy gap between the lowest doublet of states of Ho{S), 
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and the higher excited states. Thus the result ( |l.2| ) can start to fail if either the mean hyperfine coupling Wq, or the 
number N of nuclear spins, becomes too large. 

(iii) The "bias" energy = Je{Si — S2) ■ Ho, where Si and S2 are the two semiclassical minima of Ho{S), must 
also satisfy £,h -^^o- 

Conditions (ii) and (iii) simply prevent the two lowest levels of Ho{S) from significantly mixing with higher levels, 
once the nuclear spins or magnetic field are coupled in. 

Although these results are derived for the particular example (1.3), we believe that the instanton method we use 
can be applied to most reasonable forms of H{S, {(?}) like ( |l . l| ) , to give a truncated form like (1.2), provided the same 
"low-energy" restrictions (i)-(iii) are satisfied; all that will change will be the values of the parameters in (1.2) (for 
the case (1.1), the values of the parameters appearing in ( |1.2| ) are given in Sec. Ill below). 

In what follows we begin (in section ||) by describing the central spin model, and giving some discussion of the 
relevant physical coupling energies in it. In section [II we show how an initial microscopic Hamiltonian, describing the 
giant central spin and its microscopic environmental satellites, can be truncated to a low-energy effective Hamiltonian 
description. In section IV we proceed to the discussion of the results of the exact diagonalization studies and, finally, 
in section we summarize our results. 



II. THE MODEL, AND SOME SIMPLE EXAMPLES 

In this section we b riefl y describe the physical origin of the central spin model for nanomagnetic systems, and the 
Hamiltonian of eqtn. ( |l . l| ) , and also remark on what is left out of it. Since the variety of physical examples is rather 
large, as is the range of possible coupling energies, we also mention some reference examples. 



Let us consider first the central spin itself, whose dynamics is described in eqtn. (LI) by Ho{S). A typical central 
spin is made up of a very large number of microscopic spins, whose motion is assumed to be locked together by 
nearest-neighbour exchange, or "Hund's rule", couplings (or, in the case of superparamagnetic spin complexes, by 
indirect spin-spin couplings such as the RKKY interaction). These couplings can be very large, as much as 0.1 — 1 eV; 
this is much larger than other individual spin energy scales such as those coming from anisotropy, so that in principle 
one can have large monodomain magnetic particles, with many spins (up to perhaps 10* , depending on the particular 
systems) lined up to form a "giant spin" p|J^,p^. Another more subtle example is the antiferromagnetic grain, in 
which nearest neighbours are anti-parallel, and one obtains a "giant Neel vector" |lj,|l^. The superparamagnetic 
spin complexes are more delicate still, and typically arise at low temperatures when paramagnetic electronic spins in 
a disordered magnet "lock together" via indirect exchange, RKKY interactions, etc 

The first step in analysing such a system is to treat the entire spin complex as a rigid quantum rotator 5, with 
dynamics governed by an anisotropy field in the particle/grain; i.e., write a central spin "bare Hamiltonian" Ho{S), 
with I S* 1= 5 a constant. The next step to take, if the temperature is low, is to try truncating a Hamiltonian like 



(1.3) to a 2-level system, involving only the lowest 2 levels. This has been done by various authors; in fact the papers 



of van Hemmen and Suto ^ analyzed a much more general class of bare central spin Hamiltonians than (1.3). 

The giant spin model omits the higher excited states of the particle, involving spin flips of spins inside the grain 
(internal magnons); the neglect of these is usually justified by appeal to the large value | J | of spin exchange P,p^, 



in models like equation (1.3). Incorporation of such processes implies relaxation of the constraint | S Inconstant. In 
most circumstances this will hardly be a serious omission; the lowest "internal magnon" excitations do not even reach 
the top of the "ginat spin" multiplet of (2s -I- 1) levels until the grain or molecule reaches a radius L ^ (RaY^^, for an 
AFM grain {R is the domain wall length, roughly equal to the maximum size of a monodomain magnet; and a is the 
typical interspin spacing in the magnet); or L ^ {R^a^Y/^ for a FM grain. Even when we do get mixing of internal 
magnon modes and giant spin levels at these energies (i.e., energies ~ 2SK, where if is a typical anisotropy term), 
this will hardly affect either the energetics or dynamics at energy scales flo ^ K or below. Another omission is that 
of surface excitations of the system - in any real nanomagnetic grain, there can be low-energy magnetic excitations 
involving "loose spins", i.e., spins which, because of imperfect preparation or the inevitable defects, dislocations, etc., 
have a coupling energy ^ J to their neighbours. These couplings are hard to quantify, but it is rather unlikely that 



they will be as small as the single- ion anisotropy energies in (1.3). We will also ignore the possibility that surface 
spins might have anisotropy energies differing from those inside the grain. 

Finally, in this paper we will ignore all couplings to external bosonic baths such as electrons, phonon, or photons. 
Electron effects on nanomagnets were recently discussed by us |l6|] ( as well as by Fukuyama and Tatara for domain 



walls [^). Phonon effects have been discussed by many authors ||2l|,|2^. Photon effects were shown to be negligible 
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some time ago Experimental predictions for the behavior of nanomagnets coupled simultaneously to oscillator 
and spin bath were given in Ref. ; the work was also recently applied to the Mni2-acetate system |Q . 

We now turn to the problem of main interest here, i.e., the couphng of S to other spin degrees of freedom - as 
already noted, we are primarily interested in nuclear spins, both inside and outside the grain. A nucleus with a finite 
spin /, and spin moment —jn^I, interacts both locally with the electron clouds at the same ionic site, via the contact 
hyperfine interaction, and also non-locally with other ions via the dipolar field. There may also be other residual 
interactions such as transfer hyperfine interaction (as for, e.g., F^^ in MnF2 or C0F2), or quadrupolar coupling if 
I > 1/2. Finally, one may have interactions between the different nuclear spins, via, e.g., spin waves (the Suhl- 
Nakamura interaction), or dipolar couplings - these are very weak (Suhl-Nakamura interactions are ~ 10~^ and 
nuclear dipole-dipole interactions ^ 10"^ K) but, it turns out, can play a role in the dynamics of S. 

The most important interactions are the contact hyperfine interaction —"fffhl ■ A - So to the electronic spin on the 
same ionic site, and the dipolar interaction between the nuclei and electrons. The strength of these interactions is 
measured directly in NMR experiments; the hyperfine couplings range from 1 — 3 MHz (i.e., 5 — 15 x 10^^ K) for 
protons in H^, up to values greater than 5000MHz (0.25 for some rare-earth magnetic nuclei (e.g., Tb^^^, Dy^^^); 
these latter correspond to local fields acting on the nuclei as high as 500 Tesla, and come overwhelmingly from the 
contact interaction. Thus when a central spin rotates, with all internal electronic spins rotating in unison, the contact 
interaction tries to force the nuclei on the magnetic sites to follow. If the contact hyperfine interaction strength is 
ujo, then it is clear that the ratio ujo/ilo, with fl^^ the timescale of central spin rotation, is going to be crucial to the 
nuclear spin dynamics. The parameter Qq for the ferromagnetic grain is defined by the anisotropy energy K, 
while for the antiferromagnctic grain it is much higher, typically Q ^ 2\/ KJ [|l^,|l^. If the electronic spin frequencies 
flo ^ LUo, then the nuclear spins will experience a "sudden" perturbation, and few of them will follow the central 
spin; conversely, if ujq ^ fio, the nuclear spins will follow adiabatically. In fact, the ratio ujo/^o is usually somewhere 
between 0.001 — 0.1. On the other hand for nuclei in non-magnetic ions, luq/^o may be much smaller, and will of 
course depend strongly on the host (in a magnetic host, transfer hyperfine interactions —"fN^I • Aqj ■ Sj from nearby 
magnetic moments Sj, can greatly increase cuo', for example, uJo for F^^ in MnF2 is 160 MHz). 

A large central spin will also interact with nuclear spins in the surrounding medium (a substrate or solvent) via the 
dipolar interaction generated by the central spin dipolar field. This interaction may not be negligible; for a central 
spin with S = 10**, the dipolar coupling to a nucleus at a distance of lOOA is already ^ 1 MHz, rising to 30 MHz at a 
distance of 30A; and there is clearly a very large number of nuclear spins within lOOA of the central spin! We should 
also note, that the same considerations may be also applied to quite distant paramagnetic impurities. For S — 10'^ 
they couple to the giant spin in the MHz range up to distances lOOOA. 

The approach we shall take is to include first the overwhelmingly dominant hyperfine interaction in the coupling 
between S and each nuclear spin. Since the individual electronic moments are locked together this leads directly to the 
second term in (^]^). We will drop the internuclear couplings; these are so weak that their effects can be handled in a 
very simple way once the main "central spin" problem has been solved, as a weak after-effect p^ ]. In fact it is trivial 
to show that the effects of the internuclear dipolar coupling can be entirely handled, in this effective Hamiltonian 
approach, by adding back the dipolar coupling to the low-energy effective Hamiltonian in exactly the same form as 
it appears in the high-energy Hamiltonian- this is what we ha ve d one in writing equations (1.1) and (1.2). Thus, we 



shall use, as our starting point for analysis, the model in Eq. (1.1), without the internuclear dipolar interaction 



III. THE EFFECTIVE HAMILTONIAN 



The task of this section is to truncate the "giant spin" description, and find a new effective Hamiltonian which 
operates only at low energies. In fact we reduce the (25* -I- l)-dimensional Hilbert space of the giant spin itself down 
to two lowest levels - the giant spin is now described by a Pauli vector matrix r. This two-level system interacts 
with a "spin-bath", described by Pauli matrices {(Jk}, with k ~ 1,2,3, .. .TV. The new effective Hamiltonian thus 
operates in a reduced Hilbert space, of complex dimension 2^+^. We may use this Hamiltonian H^ff provided higher 
energy degrees of freedom are not excited. This means in practice that H^/f is restricted to energy scales ^ fto, the 
characteristic frequency of zero point fluctuations of S. The higher levels do of course affect H^ff', in the usual way, 
their effects are incorporated into renormalised "coupling constants" in H^ff. 

In what follows we first explain how the truncation procedure is done for this kind of problem; then we derive H^^ff, 
starting from the giant spin model. 
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A. Free Giant Spin 



The technique we use is the standard one of separating the " slow" and " fast" degrees of freedom in the problem; for 
our giant spin the difference between the two is illustrated in Fig. 1. The truncated Hamiltonian should only include 
slow variables, the fast ones having been eliminated in favour of a new set of couplings between the slow ones. 

Let us first briefly recall how the truncation is done in the absence of the environmental spins and magnetic field - 
this makes the method clear, and also gives us some results we will need later. Since 5^1, one uses semiclassical 



methods - the task is to start from a free giant spin Hamiltonian Ho{S) like (1.2), and derive a reduced Hamiltonian 
H°j:j{t), operating in the space of the lowest two levels of Ho{S). We shall choose a basis in this truncated space 
such that the eigenstates of correspond to the 2 semiclassical minimum states of Ho{S), defined by coherent state 
vectors |ni) and |n2), such that {ni\S\ni) = Sfti and {n2\S\n2) = Sn2; the eigenstates of H°jj(f) are then linear 
combinations of and 1^2), which we can determine once we have found the four matrix elements {na\H°j j\fii}) 
with a, (3^ 1,2. 

Formally one can do this as follows, for the free spin. Consider the path integral expression for the transition 
amplitude r°^(i), during the time t; this is given by [5[p^: 

rn(T=t)=n^ f ft > 

Kp{t) ^ ^ Vn{T)e^p\- dTCo{T)\ , (3.1) 

Jn{T=0)=n^ {Jo J 

where the free spin Euclidean Lagrangian is 

Co = ~iS0ipsin0 + Ho{n) . (3.2) 

Here 6 and cp are the usual polar and azimuthal angles for the unit vector n^r). 

Now in the semiclassical approximation there are two fundamental time scales in the paths n{T) in ; these are 
r2~^, the time required for the instanton traversal to be made between states, and A~^, the typical time elapsing 
between instantons. By definition, an effective Hamiltonian is supposed to reproduce the slow dynamics of the system 
in the truncated Hilbert space of the two lowest levels, i.e., for long time scales an evolution operator is approximated 
as 

Kp(.t) « (e-^^="*)„^ . (3.3) 
Since Aq is exponentially smaller than flo, and nondiagonal elements are i?e// ~ Aq, we can write 

« S^p - it{no,\H:^j\n0) ; {^-^ « t « A'') ; (3.4) 
Then we immediately find the matrix elements of iJ°^j(T) for a ^ f3 as 

{Kff{T))^p = \Kp{t) ; (f7-i « t « A-i) . (3.5) 
As a concrete example, consider the easy-axis/easy-plane Hamiltonian ( |1.3| ), where 



sin^ 9 + sin^ 6 sin^ if 



(3.6) 



The two lowest states are fii = z and fi2 = —z; henceforth we write these states as | ff) and | 1J-). In the usual case 
where K±/K^^ ^ 1 (so that the tunneling amplitude is appreciable) one has only small oscillations of (p about the 
semiclassical trajectories tp — or tt, and by eliminating ip one has 

Coie) = -^0' + SK\\ sin^ e , (3.7) 

giving a classical equation of motion = Clo sin 6', and instanton solution ||^,|l^, going from | fl") to | 4) 

sin6'(T) = l/cosh(r2oT) (3.8) 
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(centered at r = 0); in this system 



1/2 



(3.9) 



The frequency flo then sets the ultraviolet cut-ofF for the Hilbert space of ff°jj(r), and one finds, by substituting the 
semiclassical solution into (3.1) and evaluating a determinant over the quadratic fluctuations around the scmiclassical 
solution psl (the zero mode contribution in the usual way gives a factor it), that from (|3.5|) we get pljlF 



(3.10) 



^o{S) = -Y^ \J -ReA^o'^VLoeyi^{~A^^^ = 2AoCosttS , 



(3.11) 



(3.12) 



where the action aI^'' is that for transition between the two semiclassical minima, cither clockwise ( 77 = +) or 
anticlockwise {rj = — ); the phase rjnS is the Kramers/Haldane phase, coming from the linear in time derivatives 
kinetic term in (3.2). Without this phase, we would simply have a splitting |Ao| = y^2ReAo/'!Tfto exp{— Aq} with 
Ao - 2S{K\\/Kj_W2_ 



B. Giant Spin in External Magnetic Field 



Let us now include a weak magnetic field {h = je \ Ho \ /2i^|| <C 1). The effect of an applied Ho is mostly nicely 
seen in the instanton language; since the 2 possible paths between the degenerate minima involve opposite Haldane 
topological phase |2^,0 , this phase can be changed by an external field j2^ , causing the tunneling splitting to change 
and even oscillate. Then the corresponding Lagrangian written in terms of angles {9, if) has the form (in imaginary 
time) 



C ^ Co + 5Ch 
5Ch = -'jeS{H^ sin 9cos(p + H^ sin 9sm(p + cos ( 



(3.13) 

where Co was defined in ( |3.2D . As before, we perform the Gaussian integration over the small deviations of around 
zero or tt, and obtain the effective Lagrangian for 9. Keeping only linear terms in magnetic field we write 



5Ch{9)^--i,S{H^oSm9 



2Ki 



(3.14) 



whereas Co{9) has the form (3^). With the same accuracy in /i <C 1, the correction to the effective action is given by 



where So^\t) is the semi-classical trajectory defined by (3.8). Substituting (3.8) in(3.15) we find 



Mi") = -V^{^+^^K^^/K^iJ) ■ Ho , 
and the final expression for the effective action in the weak magnetic field takes the form p£ 

7e7rS' 



A(^)=25J^ 



(3.15) 



(3.16) 



(3.17) 



If Ho is directed along z (easy-axis), this biases the symmetric 2-well problem in (1.3). Within our accuracy (we kept 
only linear terms in Ho) the correction to the topological phase due to H^ iz zero. However if applied along x, the 
magnetic field lowers the barrier height and displaces the degenerate minima towards each other in the XZ-plane. 
In this case the clockwise and anticlockwise paths become inequivalent and the tunneling splitting changes with the 
magnetic field. Application of Ho along y lowers the barrier and distorts the semi-classical path. This field leads to 
oscillations in the value of the tunneling splitting since the second term in the parentheses is nothing but the 
Berry phase accumulated by the giant spin during the tunneling. 
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C. Giant Spin plus Spin Bath 



We now include the interaction of the giant spin with the environmental spins. As discussed in section II, we start 
with the Hamiltonian ( [l.l| ) (after dropping V^ip), where the {uJk} are the hyperfine couplings at each site, and the 
N Pauli variables {(Tk} describe each of the spin-1/2 bath modes. Without the hyperfine term, the total nuclear 
spectrum, containing (2)^ lines, is almost completely degenerate, with only a tiny spreading ~ Tj"^ of levels caused 
by the internuclear dipolar interaction Vdip- With the hyperfine interaction, the nuclear levels are now spread over a 
total range ~ ujoN (assuming ujk ~ for all nuclei), with most levels concentrated in a Gaussian peak of half- width 
- Since cuo » T^^, this means that the nuclear spin spectrum and dynamics are basically slaved to the 

spectrum and dynamics of S. This of course sharply contrasts with the usual oscillator bath environments, where the 
bath spectrum is only very weakly perturbed by the coupling to a macroscopic system. 

What we wish to do is to truncate to a low energy effective Hamiltonian _ff°yj(r, {<Jk}), where r, as before, operates 

in the subspace of the lowest two levels of S. The general way to do this is to separate "slow" (frequency ^ flo) from 
fast (frequency S> ^o) processes, in the combined system/environment dynamics. One then incorporates the fast 
processes into parameters of the low-energy effective Hamiltonian. In the instanton method, there are two kinds of 
processes that one may consider, shown schematically in Fig. 2. In the periods between instantons (i.e., when n — fii 
or n2), there are diagonal couplings between the environment and the system - in general, if one takes into account 
the bath dynamics, these will be dynamical. On the other hand, the transition between ni and n2 (the instanton) 
can obviously perturb the environment - thus there must, in the reduced Hamiltonian, be a dynamic non-diagonal 
coupling generated between the system and environment. The large values of cok mean that during the instanton 
transitions of 5, a great deal can happen in the nuclear bath - the time- varying hyperfine field LUkS(T)/S, acting on 
each (Tk during the transition, can cause at to flip. Since one typically has ujk/^o ^ 1 (weak couphng regime), the 
probability that a particular dk will flip during a single instanton passage is\ak\'^/2, where \ak\ '^-^ najk/^flo (this 
result is derived below, in the course of our derivation of Heff, see also refs. |0,§|). Thus in each tunneling transition, 
roughly A = 1/2 J2k I'^feP nuclear spins flip, and often A > 1. Thus in the couphng to the spin bath we must deal with 
dynamic interactions between S and the dk , during the instanton itself, when making the truncation to Hgf f . These 
will yield terms non-diagonal in the giant spin basis, in which r+ or t_ may be coupled to several environmental spins 
at once. 

Diagonal coupling: We first deal with the truncated static interaction which must exist between r and the 
{dk} when no tunneling is taking place; this is of necessity a diagonal term. Let the two relevant orientations of S 
are Si and 5*2, with 5*1 = —S2- Then the only interaction with the spin bath is the hyperfine interaction, and the 
static contribution to the effective Hamiltonian is 



Heff^lr^T^^kdl, (3.18) 
fe 



in the reduced Hilbert space. 

In the more general case where 5*1 and S2 are not antiparallel, we can still define a basis where |ni) — \ fl") and 
1^2) — I -D-)- ^ow for the two different orientations Si and S2, the effective fields acting on dk will be 7^^'' and 
These fields are the sums of the hyperfine fields L0kni,2 and any other static fields acting on the nuclei. We define the 
sum and the difference terms as 

II r -.(1) -.(2) 
^kl-k ^Ik -Ik 

u^irhu^-f^^^+^i^K (3.19) 

where the Ik and are unit vectors. Then the truncated diagonal interaction takes the form (we project on states 
I fl") and I JJ.) using standard (1 + fz)/2 and (1 — Tz)/2 operators) 

jv, 1-^ I ( ^ ^ ^ 

H?ff - E -^+t'-^]-'-^^\\-^Y: -I • + E ^4 ' (3-20) 

k=l ^ ^ ^ k=l k=l ^ 

i.e., one term which changes when S^i — > S2, and one which does not. Usually lu^^ ^ lo^ , and lu^I ^ lUq. 

Non-diagonal coupling: Let us now turn to the dynamic interactions which occur when S is tunneling. There 
are two effects we must deal with - the effect of the nuclear spins on the giant spin dynamics during tunneling, and 
also the effect of the motion of S on the nuclear spins themselves. In general both effects have to be handled together, 
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in a mutually consistent way. In this paper we concentrate on the weak coupling regime U)k/^o *C 1 when the idea 
of an effective Hamiltonian is meaningful (we recall that by definition H^f / is operating in the low-energy subspace, 
and that a consistent solution requires ujk ^ f^o)- It is then fairly straightforward to handle interaction effects, in an 
expansion in Uk/^o- 

Since we want an effective Hamiltonian which describes both f and the {<Jk}, we examine, instead of the transition 
amplitude in (^), the more general matrix element 



N 



k=l " 



"3 



'Dn{T)exp ■ 



I dT[£<,(r)+^a(r)+<5/:.(r)] 

k=i 



(3.21) 



'Cfc(T) = --6'fc.^fesin6'fc 



(3.22) 



N 



5C„{t) = ^CTfe(r) • 7fe(r) 
where 7fe(T) is now the time-dependent field 



fc=i 



Iki^) 



^kSjr) 
2S 



(3.23) 



(3.24) 



and Co{t) was defined in ( |3.2[ ). We assume in this calculation that the instanton is fast and the nuclear dynamics 
is slow and therefore we first solve for the instanton trajectory minimising Co{t) + SC{t) and calculate then the 
transition amplitude from \na > to Inp >. Since we deal here with non-diagonal matrix elements (a ^ /3), we deal 
only with the actual transitions, i.e., with the matrix elements T^j,-^. 

The result of such a calculation will depend, of course, on the particular form of the giant spin Hamiltonian, the 
spin-bath Hamiltonian and the interaction potential. Let us therefore perform the calculation for the easy-axis/easy- 
plane model (l.S) in zero external magnetic field coupled to the nuclear spins via hyperfine fields. The coupling term 
( 3.23 ) in the Lagrangian then takes the form 



N 
k=l 



(r) sin 9 cos ip + a'^ (t) sin 8 sin ip + (r) cos 6 



(3.25) 



The analogy between this term and (3.13) is obvious. Within the same accuracy, the correction to the minimal action 
from the hyperfine coupling is given by 



N 

2 

fc=i 



fc=i 



where the instanton-generated field acting on the nuclear spins is defined as 



h^^^fir) = (?7sin6'('')(T), zyy^xi/^sin 6'('') (r), cos6'('')(t) 



(3.26) 



(3.27) 



For the microscopic nuclear spin 1/2 we do not use a quasiclassical method to evaluate the transition amplitude. 
Instead, we notice that the problem of evaluating the dynamics of X]fe['^fc(''') + ^^'k'^i'^)] identical to solving the 



spin rotation of each of the nuclear spins in the time dependent magnetic field ujkh^f f{T) /2. The integrals ( 3.26 ) are 
easy to evaluate for the instanton trajectory in the easy-axis/easy-plane model and we find [compare with ( 3.17 )] 



N 

k=l ° 



(3.28) 



A note is in order here. We do not use the Wick rotation from imaginary to real time to calculate the evolution 
of the microscopic spin wave- function in the time-dependent instanton field 7fc(T). To justify this we consider the 
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text book solution for a related problem. It is formulated as follows. A spin-1/2 particle passes through a 1-D 
barrier along the y - direction. An external magnetic field is nonzero only inside the barrier. What is the spin 
structure of the transmission coefficient (say in the basis set of eigenfunctions of CTz)? It is easily shown that one has 
T = To exp{— ckfTa;}, where To is the transmission coefficient when the external field is zero, and a (which depends on 
Hx) is real. For this reason, we conclude that the co-flip amplitudes have to be calculated in imaginary time, which 
means that the evolution of the environmental spin does depend on whether the external field acting on it is classical 
or originates from the tunneling motion of a combined system. 

To get the final form of the non-diagonal terms in the effective Hamiltonian we sum up the contributions from 
the external magnetic field and from the hyperfine coupling to the spin bath and write the result for the transition 
amplitude as 



f m(0 =^tY. \[^R^o^o exp I -4") - ir^Y, 

r]=± I fc 



akU- ak + irif3on ■ Ho 



(17„-i«t«A„-i); 



(3.29) 



Qffcn : 



IX, 



K\\/K^ y) ■ 



(3.30) 



l3on = -^{-ix, ^K\i/K^y), 



(3.31) 



where T^j^-ff is an operator in the nuclear spin subspace. We may now give the result for the non-diagonal coupling 
term in this easy-axis/easy-plane model 



ND 
eff 



t 



T^ti^^it) + H.c. 



T- Ao ^ exp |-i?7 yrS + ^ akfi ■ at - (ion ■ Hoj j + H.c. 



= 2A„f_ cos 



vrS* + ^ akfi ■ dk - Pofi ■ Ho 



H.c. , 



(3.32) 



where Ao ~ —^j2ReAol'K i7o e ^'^^o as before. Adding the diagonal term (3.20) we get the final truncated effective 
Hamiltonian: 



H^ft = < 2AoT_ cos 



7r5 -f ^ aun ■ dk - l3on ■ Ho 



fc=i 



H.c. 



, N ^ N 



fc=i 



(3.33) 



with the parameters Ik — z, lu^^ = ujk, ^k ^ ^'^^ easy axis/easy plain Hamiltonian ( |l.3D ; the parameters ak, ftk, 
and (3o are given in (3.30) and (3.31). 

This result for the easy- axis/easy-plane model is not the most general result we can get for the problem of a giant 
spin coupled to a nuclear spin bath. More generally we will find that both the tunneling splitting Ao and the phase 
T^S will be renormalised by the coupling to the nuclear spins - these renormalisations are 0(w|/rio) for each nuclear 
spin (they correspond to a change in the instanton trajectory caused by the couphng to the {tTfc}, and hence appear 
first in the 2-nd order in this coupling). Thus more generally we find that ttS — > $, and Ao — > Ao, where 



(3.34) 



Ao = Aoexp{-^4} 



(3.35) 
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Physically, the {4>k} are the Berry phase terms to be added to the instanton phase, coming from the nuclei; and the 
{6k} describe the renormalisation of the tunneling splitting. We do not attempt to calculate either ipk or Sk here, 
since their effects are not going to be observable - we only observe Aq and not Aq or the "bare" phase nS. 

We have now achieved our main aim in this section, of reducing the giant spin Hamiltonian to a truncated one. It 
is interesting to notice what kinds of process are now included in this effective Hamiltonian. Notice that if we expand 



out the cosines in (1.2), we see that we have a whole series of terms like 

T±r„^^5...<<a^3<-- - . (3.36) 

in which the instanton flip of the giant spin couples simultaneously to many different nuclear spins. Thus, as well as 
the simple diagonal and non-diagonal transitions, we also have the possibility (see Fig. 2) of multiple transitions in 
the bath, stimulated by a single instanton. The parameter A = ^ la^p - as already noted, measures the average 
number of nuclear spins to be flipped during each instanton. 



IV. EXACT DIAGONALIZATION STUDIES. 



In this section we wish to verify our analytical derivation (|1.2| ) of i?e// by an exact diagonalization method, and 
to analyse the structure of the low-energy effective Hamiltonian for arbitrary values of parameters S and K±/K\\, 
including regimes where the WKB approach breaks down. 

For simplicity we consider the Hamiltonian for a giant spin coupled to a single nuclear spin 1/2 (for small ujk/^o 
nuclear spin contributions are additive); a similar tactic was employed by Abarenkova and Angles d'Auriac |p^ . We 
have 



Si + AS, 



-S-ff 



HqS 



(4.1) 



where we measure energy in units of K\\^ magnetic field in units K\\/^e, and we define the anisotropy parameter 
A = K^/K\\. Let us also apply the magnetic field Ho along the y axis (i.e.. Ho — Hy). For this particular choice of 
the external parameters we can rewrite (1.2) as follows 



H 



— < 2Aof_ cos 



7r5 — ip + an ■ a 



H.c. '> + 



LOoHy 



(4.2) 



where the last two terms describe the static diagonal interaction (and as before, we neglect both the interaction 
between the nuclear spin and magnetic field, and the internuclear dipolar coupling). In this equation we have defined 



an = ^ ( - ix, y/ VA) , 



(4.3) 



nSHy 

lb = — ^ 

^ 2A 



(4.4) 



Let us start by defining the diagonalization procedure. Since the giant spin only weakly couples to the nuclear 
spin 1/2, we hav e four lowest energy levels, well-separated from higher levels by the instanton frequency fio- The 
Hamiltonian (42) is given in the Hilbert space of these four levels. Therefore, we can write 



{a\H,ff\b)^{a\H\b); (a , 6 = 1 , ... , 4) 



(4.5) 



or 



(a|Fe//|&) « (a|5]£;^|*^)(vl/^|6) = 5^(?7+),^£;^?7^, , 

7—1 7—1 



(4.6) 



where 



(4.7) 
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and the |^^) are the four eigenfunctions of the Hamiltonian (4.1), corresponding to the four lowest energy levels 
El , ... , Ei. In fact, Eq. ( |4.5| ) is nothing but the definition of the effective Hamiltonian. 

The choice of the appropriate basis set of \a) is important here. The most convenient and physically transparent 
basis set is that of the nonintcracting system, i.e. the eigenfunctions oi H = Ho{S). Let 1$°) and |$2) be the 
eigenfunctions corresponding to the two lowest levels of Ho for the giant spin 5*. Then the functions 



\X±) 



V2 



(4.8) 



are localised single- well states, which in the semiclassical limit correspond to the previously introduced \ni) and |r?2). 
With nuclear spins added we define the basis set \a = 1, 2, 3, 4) as 



|1> = 1X4 



|2) 



\X+) 



|3> 



IX- 



, |4> 



IX- 



(4.9) 



(where (g) stands for the direct product of state vectors). This particular choice of the basis set \a) is the starting 
point for our diagonalization procedure. 

The diagonalization procedure itself is arranged as follows. First, one calculates the eigenstates of the nonintcracting 
problem to construct the set |a), and then diagonalizes the interacting Hamiltonian to obtain the energy spectrum 
and the unitary transformation U as defined in Eqs. (4.6) and (4.7). As H^ff is a matrix (4 x 4) matrix it can be 
written as 



Heff = 



3 



'3 ' 



where fi and 2 ■ ai (i = 1, 2, 3) are the Pauli matrices and Tq, cTq are the unit matrixes. Thus 



Tr{H,ff{n(g>aj)} 



(4.10) 



(4.11) 



We analyse these coefficients and compare them with those resulting from (4.2). If S is an integer we get (we omit 
here the (8) sign) 



(—1) < 2A0COS?/; • "TrCTo — 2Aoasmip ■ fy&x + 2Ao— = sin?/; • fj.(Ty 

vA 



LUoHy 

4A 



TnCTy 



TTUJo 
2K 



(4.12) 



and for S = 1 /2-integer we get 



Heff « (—1)^'^'^"'^ '|2AoSin'0 • Tx&o + 2AoQ;cosV' ■ fyO-x — 2Ao 



VA 



COSTp ■ TxOy 



(4.13) 



Now we turn to the results. 



(a) Semiclassical regime S 3> 1: We start by carrying out the numerical analysis for a wide region of parameters S 
and A and as an example we present our results for S = 50, A = 50, and lJo = 0.2 (recall again that energy is measured 
in units of K\\). In full agreement with (4.12) the effective Hamiltonian H^jj, derived by the exact diagonalization, 
has only five non-zero coefficients C^, i.e.. 



H. 



ED 



Cxo ' '^x^o ~(~ Cxy ' '^x^y H~" Cyx ' '^yO'x ~h Coy ' ^o^y ~\~ Czz ' '^z^z 



(4.14) 



In Fig. 4 we present the behavior of the 4 non-diagonal coefficients as functions of the magnetic field Hy, and compare 
with the corresponding analytical forms (4.12), coming from the instanton analysis. As to the value of Cz^i numerically 
it is 0.0972 in comparison with the analytical result 0.1, showing a small zero point vibration (z.p.v) reduction of 

(x+|^.lx+) < s. 

We also have performed calculations for half-integer spins. For S = 49.5, A = 49.5, and lOo — 0.2, for example, we 
have found a similarly good agreement between the numerical and analytical (4.12) results. 
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Bearing in mind that we kept only the leading terms in the hyperfine coupling when deriving (4.12, 4.13) we 
conclude from these results that the effective Hamiltonian (|4.2|) (and (1.2)) gives a complete description of the low- 



energy dynamics of the system in weak magnetic field for S ^ 1, and that moreover the values of the coefficients 
computed for the effective Hamiltonian by the instanton expansion are in rather good agreement with the exact 
results. 

(b) Beyond WKB: Arbitrary values of S and A: Although the above results strongly support the validity of 
our effective Hamiltonian (^]^) for S* ^ 1, it is of considerable interest to explore its validity when we are far from 
this WKB limit. In fact it was shown in the original papers of van Hemmen et al., as well as Enz and Schilling 
that WKB results for the free giant spin are remarkably accurate even when S is less than 10, at least as far as the 
calculation of tunneling splittings are concerned. Thus we shall now use our exact diagonalisation results to tell us 
how accurate is our H^ff in (1.2), when S and A are arbitrary, and the action Ao is not much larger than unity. 



We begin by defining a function "^{S, A, Hy), such that the actual tunneling splitting Ao(^') = AqCOS^' for the 



ip; we wish to see deviations for small S and 

like cos^ (and 



non-interacting problem (for 5- integer). In the semiclassical limit 4* 

large A. Notice that even for small values of S, the coefficient Cxo in (4.14) still oscillates with Hy 
Cxy and Cyx like ztsin^*). 

In Fig. 5 we present our numerical results for ^E"/?/; as a function of S and A, with H" equal to some value smaller 
than one period of oscillation of Ao(\l/) [in fact, in most cases ^ ^ Hy rather precisely unless 5* < 3]. Wc plot the 
ratio \E'/V' as a function of S at A = 50 and A = 10 (Fig. 5. a) and as a function of A at 5* = 50 (Fig. 5.b) to see where 
the analytical prediction (4.4) becomes invalid. As we can see from Fig. (5. a), the difference between the exact value 
of 5* and the semiclassical result 'i' — tl; occurs at small values of S, whereas for S* > 10 this difference has almost 
disappeared. Moreover, the curve at A = 10 is falling to 1 more rapidly than at A = 50, because with increasing 
effective action Ao the WKB description becomes more accurate. Fig.(5.b) shows the violation of the WKB regime 
with increase of A at 5 = 50 as the effective action Ao goes down. At S* = 50 the analytical formula (4.4) for i/j begins 
to show appreciable deviations from the exact result once A > 700. 

So far so good; we now turn to the analysis of the interaction Hamiltonian itself, by considering the dependence of 
the coefficient a on 5 and A. First, we calculate a as a function of S* at a constant value of the effective action Aq. 
One expects that for large Ao there should be good agreement between the numerical result and the analytical result 
(4.3). In Fig. 6 we plot the coefficient Cyx, divided by 2AoSin^', together with a, defined in (^^), vs the giant spin 
quantum number 5* at A = 5*^/50. With this definition of A the effective action Ao is a constant « 14. Since Aq and 
^'(S', A, _ff°) are already known numerically, we plot the results of the exact diagonalization for a together with the 
analytical one. From Fig. 6 we see that once S > 10, the result (O) for the effective coupling a is extremely good. 

To study the deviations from (O), we show in Fig. 7 numerical results for Cyx / {'2 Ao sin together with a from 
(4.3) vs the giant spin quantum number at A = 50 and A = 10 (Fig. 7(a)) and vs A at 5' = 50 and S' = 10 (Fig. 
7(b)). The solid lines in Fig. 7(a) (parallel to the ^-axis) are the analytical results. We can see, that at small 
values of the giant spin (at 5 < 10 in this particular case) it becomes useful to derive numerical results to describe 
its interaction with the nuclear spin. Notice (Fig. 7(b)) that the numerical (triangles) and analytical (solid) curves 
coincide extremely well for S' = 50 whereas for 5' = 10 the differences between the numerical and analytical curves 
become obvious once A > 100. 



Thus we can conclude that (4.3) becomes more and more accurate with the increase of S at any fixed value of A 
once the effective action Ao = 25'/A^/^ becomes significantly greater than unity (compare, e.g., in Fig. 7(b), the point 
A = 500 for 5 = 10 (triangles down), where the effective action Ao < 1 , with the point having the same value of A, 
but S = 50, for which Ao > 1. 

Finally we analyse the effect of z.p.v. of the giant spin on the value of the diagonal coupling. The orientation 
of the giant spin is fluctuating around the classical energy minima, resulting in a reduction of the static diagonal 
interaction. It is easy to see from (3.7) that the value of < sin^ 9 l/Ao- From this it follows that the coefficient 
ujo in the longitudinal diagonal interaction (ie., that proportional to • ct^ in (4^)), can be corrected by replacing 

LUo ^ Wo 

(1- l/Ao)i/2. Thus, with increasing A, the giant spin becomes more and more " delocalized" , and the extent 
of the delocalisation is parametrized in the exact diagonalisation results by the z.p.v. reduction of Czz- To show 
how accurate is the simple replacement lUo — *■ Wo(l — l/^o)^^^, we compare in Fig. 8 the coefficient 2Czz/'^o against 
(1 — l/Ao)^^^, both plotted as functions of at A = 50 (triangles up) and at A = 10 (triangles down). 



(c) Zero magnetic field: Nonlinear corrections in a: 



Our instanton derivation in section HI was carried 

integer 



out only as far as Ist-order in the effective coupling ak- However we see immediately that in (4.2), when S 
and in zero external field, the nondiagonal term is affected by the coupling only in the second order in a, i.e., 
Aocosh(a) « Ao(l 4- a^/2). As discussed above [see Eq. ( 3.35| )] in a consistent calculation one may also have 
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corrections ~ o? which are not accounted for in Eq. (4.2). To find these corrections we have performed the necessary 
numerical calculations for for integer and half-integer S in zero external field. 

In Fig. 9(a) we show the behavior of the coefficient Cxo/'^^o (or tunneling splitting renormalisation) vs a for integer 
S. Our analysis shows that the behavior of Cxo in the region < a < 0.5 is very close to 2Ao(l - aV2)- These 
results suggest that the parameter 5 in Eq. ( |3.35| ) is given by a^, and the nonlinear correction may be accounted 
for by substituting Ao ^ Ao(a) = Ao(l — a^) in the effective Hamiltonian (p^). In Fig 9(a) we plot the numerical 
results for Cxo/'^^o (triangles up) together with the analytical approximation (1 — a^)cosh(a) (solid line). All other 
coefficients (except Czz) are zero. 

To see whether the ansatz Aq ^ Ao(a) really works well, we considered also the case of half-integer S. In Fig. 9(b) 
we plot the coefficient Cyx/'^^o vs a (triangles up) together with the analytical approximation (1 — a^) sinh(ck) (solid 
line). As we can see from this figure, our approximation works quite well in the region < a < 0.5. 

This completes our numerical analysis of the low-energy effective Hamiltonian in the problem of the giant spin 
coupled to the spin enviroment. Our main conclusion in this section can be f orm ulated as fo llows: the predictions 
about the low-energy dynamics of the system in the magnetic field, given by (O) (and (O) correspondingly), are 
rather accurate when the effective action Ao > 1. 



V. CONCLUSION. 



To conclude, it may be useful to summarize some of the esential mathematical features of our analysis. We started 



with a "high-energy" form for the central spin coupled to its spin bath environment, given by (1.1). This was truncated 



to the low-energy form (1.2) by an instanton-type procedure. However the results of this truncation were different 
from the usual quantum environment models, mainly because the spin bath modes are strongly coupled to the central 
spin. As a consequence much of the physics of the problem is contained in dynamic non-diagonal terms (in contrast 
to the oscillator bath models, where such terms are almost irrelevant), as well as in almost static but very strong 
diagonal couplings. The instanton derivation also allows one to evaluate the coefficients in the low-energy effective 
Hamiltonian, for a given initial "high-energy" Hamiltonian for a giant spin coupled to a nuclear spin bath. 

In the present paper we set out to find out how accurate this instanton result for the low-energy effective Hamiltonian 
really is. Using an exact diagonalisation procedure, we were not only able to verify that the effective Hamiltonian 
is correct, but also that the coefficients in it are rather accurately given by the "first-order" (in ojk/^o) instanton 
expansion, provided we are working in the regime where the semiclassical WKB approximation works. 

The importance of this result is simply that the low-energy effective Hamiltonian, in the form we have given, turns 
out to be remarkably easy to use in practical calculations of the spin dynamics of the central spin, as we have shown 
in a number of other papers (see, eg., ref. p6[ |, and refs therein). From the present calculations we now have a good 
idea of how trustworthy these calculations are- essentially they are as good as the semiclassical approximation itself. 
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FIGURE CAPTIONS 



Figure 1 A typical trajectory for S'z, contributing to the path integral for the propagator of a free giant spin. The 
time scale for "fast" variations of Sz while it is in one of the 2 wells or making a transition between the weells is , 
whereas the time between jumps (instantons) is much longer, roughly A~^. 

Figure 2 Processes which may occur when a 2-state system is coupled to some background environment - we show 
some possible kinds of coupling between the path of the 2-state system (solid line) and the environment. There is 
a diagonal coupling (ie., a coupling to r^), indicated by "D", to a single environmental mode (whose propagator is 
indicated by a wavy line). The first non-diagonal coupling (labelled by "ND") is to a single environmental mode, 
whereas the second one couples simultaneously to 3 different environmental modes. 
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Figure 3 Processes occuring in the propagator for a 2-state system coupled to a set of oscillator modes, after one 
averages over these modes ("integrates them out"), which allows us to join the wavy lines of the previous figure in 
the usual way. Process (a) involves the excitation of an environmental mode when the central spin is in state | 'f|->, 
and with re-absorption when it is in state | JJ.>; process (b) has the central spin in the same state for both emission 
and re-absorption. Both these processes involve only diagonal operators. Process (c) involves non-diagonal terms, 
and process (d) involves both diagonal and non-diagonal terms. 



Figure 4 The numerical coefficients Cy (4.11) of the effective Hamiltonian (4.14) in comparison with the analytical 
expressions from ( [1.12| ). The figures are: (a) Cxo (solid fine) and 2AoC0Si/' (dashed fine), (b) C^y (sofid line) and 
2Aoasinip/A^^^ (dashed line), (c) Cyx (solid line) and 2Aoasmip (dashed line), (d) Cyo (solid line) and uJoHy/AK 
(dashed line); = 50, A = 50, = 0.2. 



Figure 5 The ratio '^{S,K,H°)/^ of the exact phase ^ to the weak-field expression for ■0 (14): (a) vs the giant 
spin quantum number at A = 50 (triangles up) and at A = 10 (triangles down), (b) vs A at S" = 50 (triangles up); 
LOo = 0.1, H° = A/8S. 

Figure 6 The numerical coefficient Cyx/(2Ao sin^*) (triangles up) together with the analytical co-flip amplitude 



a from (^) (solid line) at A = S'V50; = 0.1, H° = A/8S. 



Figure 7 The deviation of a from the analytical expression (13). The figures are: (a) The numerical coeffi cien t 
Cyx / {2 Ao sin at A = 50 (triangles up) and at A = 10 (triangles down) together with the analytical a from (4.3) 
(solid lines) vs the giant spin quantum number 5, (b) The numerical coefficient Cy2,/(2Ao sin ^') at S' = 50 (triangles 
up) and at S* = 10 (triangles down) together with the analytical a (solid line) vs A; cuo = 0.1, H" — A/85. 

Figure 8 The numerical coefficient Czz/0.5 ■ uJq, vs the giant spin S* at A = 50 (triangles up) and at A = 10 
(triangles down); = 0.1, = A/8S. 

Figure 9 Zero magnetic field: (a) The numerical coefficient Cxo/^Aq vs a at 5 = 50 and A = 50 (triangles up) 
together with the analytical approximation (1 — a^) cosh(a) (solid line); (b) The numerical coefficient Cyx/2Ao vs a 
(triangles up) at S — 49.5 and A — 49.5 together with the analytical approximation (1 — a^) sinh(a) (solid line). 
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